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Abstract 

The Sivers function, an asymmetric transverse-momentum distribution of the quarks in a trans- 
versely polarized nucleon, is calculated in the MIT bag model. The bag quark wave functions 
contain both S-wme and P-wave components, and their interference leads to nonvanishing Sivers 
function in the presence of the final state interactions. We approximate these interactions through 
one-gluon exchange. An estimate of another transverse momentum dependent distribution h^ is 
also performed in the same model. 
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The measurements from the HERMES, SMC, and JLAB collaborations show a remark- 
ably large single-spin asymmetry (SSA) in the semi-inclusive processes, such as pion pro- 
duction in 7*p — >■ irX, when the proton is polarized transversely to the direction of the 
virtual photon [H Q, 0] . If the underlying process is hard, the physical interpretation of such 
single-spin asymmetry can be attributed to either the quark transversity distribution hi{x) 
01 convoluted with the Collins' fragmentation function H±{z, k±) 0, or the Sivers function 
f^{x,k±) 0] convoluted with the usual fragmentation function D{z), or both. Hi- The 
Sivers function is the asymmetric distribution of quarks in a transversely polarized proton 
which correlates the quark transverse momentum and the proton polarization vector S±. 
The nonvanishing of the Sivers function has been confirmed recently BHIiiHIil. The 
key ingredient here is that the gauge link in the gauge-invariant definition of the transverse 
momentum dependent (TMD) parton distributions f{x, kj_) generates the initial and/or final 
state interactions, which results in a phase difference in the interference between different 
helicity states of the proton 
SSA in semi-inclusive processes 
impact parameter dependent parton distribution. 

The TMD parton distribution functions are defined through the quark density matrix 

0, 




11 . There are many other approaches to understand 
15 1 . For example, in jl^ the SSA is connected to the 
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P 



-i(rfc+-l*±-fc±) 



(27r)3 

X (P5|^(r, a)^(o, ox)^(o) |P5) , 



(1) 



where is the polarization vector of the nucleon normalized to SuSj^ = — 1, p'^ is a light- 
cone vector such that p~ = 0. The gauge link C'^{^~,^±) is defined Jll . 1^ |. 



= ^exp (^-ig j v4+(?7 ,iA_)dri ^ Pexp (^ig j drj^ ■ A^{ri = oo,?7_l)^ 



(2) 



In nonsingular gauges, the second term vanishes. However, in singular gauges, such as the 
light-cone gauge, the second term will contribute. In the following, we will work in the 
covariant gauge, and so we can keep only the first term in the gauge link. 

A model calculation of the Sivers function is needed to demonstrate its existence and its 
size in the typical kinematic region. Because the Sivers function contains the interference 
between different helicity states of the proton, the model must contain 5* wave and P wave 
components to generate phase difference, i.e., involving the proton wave function component 
with nonzero orbital angular momentum 0. For example, in the quark-diquark model used 
in the proton-quark-scalar coupling contains S and P wave components, corresponding 
to the quark spin parallel and anti-parallel with the proton spin respectively. This model has 
been used to calculate the Sivers function and other interesting distributions [1, [H, 1^ 2^ - 



In this paper, we shall study the Sivers function in the MIT bag model [21[ . The bag model 



contains confine physics and incorporate SU (6) spin-flavor structure. More importantly, the 
bag model wave function has both S and P wave components. Of course, the bag model has 
a number of well-known problems, including breaking of chiral symmetry and translational 
invariance, etc.. Nonetheless, it approximately generates right the hadron spectrum (2]| : 
it yields reasonable quark distribution at low energy scale [22] ; and it can describe the 
electromagnetic form factors of the nucleon IIT 
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In the MIT bag model, the quark field has the following general forms 21 



'^aix,t)= ^ N{nK){ba{nKm)i!nKjmix,t) + dl{nKm)'ip^ri.~^jrn{x,t)} , (3) 

n>0,K=±l,m=±l/2 

where b]^ and dj^ create quark and anti-quark excitations in the bag with wave functions, 

V47r \ -a ■ xji{^i^)xm J 

For the lowest mode, we have n = 1, n = —1, and uji,-i ~ 2.04. In the above equation, 
a is the 2x2 Pauli matrix, Xm the Pauli spinor, and Rq the bag radius, x represents the 
unit vector in the x direction, and ji are spherical Bessel functions. Taking the Fourier 
transformation, we have the momentum space wave function for the lowest mode, 

I to{\k\)x m 



iV- I . 1 . (6) 



where the normalization factor N is, 

r ^ 
2Rl{uj - l)sin2cu 

The two functions ti, i = 0,1 are defined as 

1 

U{k) = j u^duji{ukRo)ji{uuj) . (7) 



It can be easily seen from the above equations that the bag model wave function Eq. (jSJ 
contains both S and P wave components. The interference between these components will 
generate a phase difference under the gauge link contribution. 

The Sivers function represents the asymmetric part of the transverse momentum 
distribution of the quark in a transversely polarized proton, and can be calculated from the 
quark density matrix A4 in Eq. (Q) through expansion 0], 

= ^ [/it(^, k±)e''''''l,P.KS, + ...], (8) 
where M is the nucleon mass. Inverting the above equation, we obtain. 



Since we work in the covariant gauge, only the first term in the gauge link Eq. Q ( the 
light-cone gauge hnk) contributes. Without the gauge link contribution, the Sivers function 
vanishes. For example, to the leading order, the above function has the form. 
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where for convenience, we have chosen a particular polarization vector Sx representing the 
proton is polarized along the x direction. Inserting the bag model wave functions Eq. Q, 
we get the following results for the leading order contribution to the Sivers function without 
the gauge link contribution, 

X y^l{hh''j^^m'{h){PSx\ai^a^m'\PSx) . (11) 
Plugging in Eq. (jSj), we find that 



tlik) + 2to{kMk)^ + tlik) 



(12) 



On the other hand, for a transversely polarized proton state, we have 

6mm' {PSx\aamO'am'\PSx) = , (13) 

as expected. 

Expanding the light-cone gauge link to next-to- leading order, we have 

oo 

x{PSx\iiAa±){^9) jdT^-A+{r.U)^l^i^am\PS.) + h.c., (14) 

«- 

where a is the fiavor index, i the color index, and T" the SUd'i) Gell-Mann matrix, g is the 
gluon coupling with quark field in the MIT bag model. 

Inserting the MIT bag model wave functions, we find that 

f^^-{xM) = -iwf^ I ^'ff)?^ ^l^A ^t + - xP^)6^^\k,^ + g^- k^){2.)5{q^) 

'Pij'Pkl {P^x I bami ^am2 ^ Pmi ^^I3mi I P^x) 



I —I 

g+ + ieq^ + ie 

[3 ,mi ,m2 ,m3 ,m4 



^vLS^lWl^Vrn2ik)ip'^m^{h)l^-i+iPm,{h - + h.C. , (15) 

where is the annihilation operator for a quark with fiavor a, helicity m, and color index 
i. In the above derivation, we have used the free gluon propagator as an approximation. 
Actually, in the bag the gauge boson propagate differently as in the vacuum M- The 
corresponding diagrams for Eq. (jl5|) are shown in Fig. 1. 
Using the identity, 

1 1 



g+ — ie g+ + ie 

we get 



z{2iT)6{q^) , (16) 



f^ri^^kx) - 2g^—j-^ J T2^'(f ^ Tip^kl{PSx\b''^miKm2^'tm-J^^Pmi\PSx) 

I3,mi,m2,m3,m4 

d^k3 ^ 
(27r)3'^™3 



X ^l^{k-q^h'^+Vm2{k) I i^,^l,{h)l'l^Vn.Ah - qi) . (17) 




FIG. 1: The leading contribution to the Sivers function fiT{x,k±) in the MIT bag model. 



The /cs integration only depends on q, and so we can write this integral as a function of g^, 
and define 

/fl^k 1 
j^s^i^^^h'^^'^rn.ih -q)^ -^F{q')6^,^, . (18) 



The function F{q^) is 

1 «, ,4 1 !• 

^r— . I Sk:, U{\k^\)U\k'.\)+ . 



T^^3o{^)i.^ - 1) 



^0(|^3|)t0(|fc^|) + ^^t0(|4|)tl(l4l) 



+ ^^l(|fc3|)to(|4|) + ^^ti(|fc3|)tl(|4l) 

where /C3 = /ca + q. It is easy to show that F{q^) — > 1 as ^ 0. 
With F{q'^), we can further simplify Eq. (fTTjl as, 



and Ci is 



Jmi+m2 



(19) 



A-!,"(x,fcx) = -V2g'C'^{S.)^ J ^F(g^)^V^L,(fc-gl)7°7V«^.(fc) , (20) 



Ca , (21) 



where 5mi+m2 factor comes from the fact that the proton is polarized along the x direction, 
and Ca is defined as 

Yl Wu^PSMn^X-rn^^l.^rn^PS.) ■ (22) 
/3,mi ,7713 

Substituting the above results into Eq. (20), we find 

Ai"(x,A:^) = -V2,U^ j ^,F{q^)~J2vl{k-q\h'l^^-U^ ■ (23) 
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From the bag model wave functions, we obtain 



V2 

k X q± ■ a 

\k'\\k\ 



idxk'Y ... ,,, (a X k)'' , ,,, 



\k' 



hik'Mk) 



\k\ 



X- 



(24) 



where k' = k — q±. Because the proton is polarized along the x direction, in the above 
equation, only the terms contribute, and the q± integral will be proportional to k^. And 
finally, we get 



/it (-^5 k_i 



TT 



[Io{x,k±) - Ii{x, k±)] 



(25) 



The two integrals Jq and Ji are defined as 



where k. 



Io{x, k±) 
Ii{x, fc_L) 
xM„ 



(27r)2 ql 
toik) + ^hik) 



k' k 

d\± F{ql)tr{k')k^-q^ 
(27r)2 ql k' kl 



(26) 



e, k = a/ {xMp — e)2 + kl with e = lo/Rq, and k' 



xMp — ey + {k± — q±y. Rq and Mp are bag radius and proton mass, respectively. In 
our calculations, we fix the dimensionless parameter R^Mp = Au. 

Working out the matrix elements of Eq. ()22|) for the valence quarks, we find 



16 



4 
9 



(27) 



for up and down quarks respectively, which means that the up quark and down quark have 
opposite signs for the Sivers function, and differ by a factor of 4. This is the result of 
the SU{6) wave function we used for the proton. For a polarized proton, the polarized up 
quark distribution has a factor of 4/3 while down quark has —1/3. Phenomenologically, 
since 7r+ production is dominated by the up quark fragmentation and tt^ is dominated by 
either the up quark or the down quark fragmentation, while vr" is dominated by the down 
quark fragmentation, the above prediction will lead to larger single spin asymmetries for 
TT"*" and 7r° than that for vr" with opposite signs if assuming the Sivers mechanism. In 
this estimate, we have neglected the "unfavored" fragmentation contribution to the pion 
production, which has been shown to play an important role for vr^ asymmetry [i^. Taking 
into account the "unfavored" {u quark) fragmentation contribution which has opposite sign 
from the "favored" {d quark) one, we will get even smaller asymmetry for vr". We note 
that the HERMES collaboration actually showed much larger asymmetries for vr^ and 7r° 
than that for tt~ [1]. On the other hand, concerning the quark distribution in the neutron, 
one shall have 4 times larger Sivers function for down quark than that for up quark by 
isospin symmetry argument from the above results, and both of them will have different 
signs compared to the proton ones. That means, with the neutron target, one would have a 
factor of 2 smaller asymmetry for 7r+ with opposite sign compared to the asymmetry with 
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FIG. 2: The Sivers functions /fip(x, k±) for the valence quarks at x = 0.3 as functions of k±, where 
the quark-gluon couphng ag = 0.2. 



the proton target. It is interesting to note that JLab will measure these asymmetries with 
the proton target, and the neutron target as well. The comparison of the SSA between the 
proton and neutron targets will provide crucial test on the Sivers mechanism for the SSA. 

In Fig. 2, we plot the Sivers functions for up and down quarks as functions of transverse 
momentum k± at x = 0.3. Here the quark-gluon coupling g is treated as a free parameter, 



and we set as = 5'^/47r = 0.2, which is smaller than the value used in 2l| to determine 
the mass splitting of baryons. Since bag model is not suitable for the calculation of the 
distribution at large transverse momentum, here we only show the results for the range of 
k± smaller than 0.7GeV. 

The contribution from the Sivers effect to the SSA in the semi-inclusive process can be 
calculated from the above results divided by the unpolarized quark distribution in the same 
model 2^1. The asymmetry Vy is calculated as 



7 X 

Vy = ^f^T{xM)lh{xM) , (28) 

where the polarization of the proton is along the y direction. We plot these asymmetries 
for the valence quarks as functions of k^^ aX x = 0.3 in Fig. 3(a), and as functions of x 
at fc^ = 0.5 in Fig. 3(b). These asymmetries are for the quark distributions. To get the 
asymmetry associated with the hadron production in semi-inclusive processes, we need to 
convolute the above results with the fragmentation functions of the hadrons. 

It is also interested to study the moments of the Sivers distribution. For example, one 
interested moment is defined as 



ft^'\x) = j d'k^ j fMx, k^) . (29) 

The numerical results for the above functions depending on x are shown in Fig. 4, which 
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FIG. 3: The bag model prediction for the asymmetry of the quark distribution in a transverse 
polarized proton as a function of and x, where as = 0.2. 
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FIG. 4: The first moments of the Sivers functions /^^ (x) for valence quarks, where = 0.2. 



can be fit with the following functional form, 

P±(i)« 



\4.06 



±{l)d 
IT ' 



x) = 0.19x^-^^(1 -x 



\4.06 



(30) 



Since the bag model is not good for small x parton distributions, we have abandoned the 
use of small x points in the fit. As an illustration, we also plot the above fit for the up-quark 
in Fig. 4. 

We can repeat the above calculations for another TMD parton distribution, hj^, which 
represents the correlation between the quark's transverse momentum and polarization in an 



8 



unpolarized proton state. It can be calculated from the expansion of the density matrix, 

M = ^[h^{x,k^)a^''k,p^ + ...] . (31) 
Inverting the above equation, we get 



hi (x, k±) = — r—— / — — -^e 

x(p|?(r,a)^^(r,a)7V75/:(o,ox)^(o)|p) , (32) 

where the proton is unpolarized. Without the gauge link contribution, this function vanishes, 
as the Sivers function does. Expanding the gauge link to next-to-leading order, we get 



2e'^k^ J (27r) 

x(P|V^.,(r,a)(^^?) / dv-Ativ'.U)ni^Y7^^aM\P) + h.c. . (33) 



where the nucleon is unpolarized. 

Using the same method as we did in the calculations of the Sivers function, we find that 

hl{x,k^) = -V2g'D^^ I ^F(g^)^¥^L(^-^L)7Vf75^™(^) , (34) 



where Da is defined as 

/3,mi?Tt3 

And finally, we can write the distribution in the form of, 

hia{x,ki_) = ^^^^^^7^^^^[Ioix,k±) - Ii{x,k^)] , (36) 

which is the same as Eq. ()25|1 except the color factor. For the valence quark distributions, 
we have 

8 4 
Du = --, D, = --, (37) 

which means that the up and down quarks have the same sign for hj^ distribution, and differ 
by a factor of two. For the unpolarized proton, the up quark distribution is two times larger 
than down quark distribution. This prediction shows that the asymmetries associated with 
hj^ for TT^ and 7r° will have the same sign. This is quite different from the asymmetries 
associated with the Sivers function we discussed before. For the neutron, one has the 
similar prediction. 

As an illustration, we plot in Fig. 5 the distributions h^{x, k±) as functions of transverse 
momentum at x = 0.3. We can also calculate the moments of the hj; functions for up quark 
and down quark, and fit with the following parameterizations, 

/i^««(x) = -1.13x1-63(1 -x)^-°% 

ht^'^^x) = -0.56x1-63(1 - x)^-°6 , (38) 
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kj (GeV) 

FIG. 5: hi{x, k±) for the valence quarks at x = 0.3 as functions of k±, where Og = 0.2. 

where the same functional dependence as f^^^ have been observed. 

In conclusion, we have calculated the Sivers function in the MIT bag model. The gauge 
link in the gauge-invariant definition of the TMD parton distribution functions plays the 
crucial role for the nonvanishing of the Sivers function. Our calculations show that the 
up quark Sivers function is 4 times larger than that of down quark with opposite signs, 
consistent with the SU{6) spin-flavor structure of the proton. These results lead to testable 
consequence for the single spin asymmetry associated with the Sivers function in the semi- 
inclusive deep inelastic pion productions: the asymmetries for vr"*" and 7r° will be larger than 
that for 7r~, and with different signs. Distribution hj^ has also been calculated in the same 
model, and we found that the up quark and down quark have the same sign, which means 
that the asymmetries associated with this distribution for vr^ and vr" will have the same 
sign. 

We end up our paper with a few comments. First, in our calculations, we have used free 
gluon propagator connecting gluon fields inside the bag, which is an approximation (2^ . 

Secondly, we have ignored the scale evolution of the Sivers function moments fyr^\x,Q'^) 
and ht^^\x,Q^) 0,123. Our calculations are performed at the bag scale, which is much 
lower tha n typic al hard scattering scales. However, the evolution of these functions is not 
clear yet [271, |28| , and is beyond the scope of the present paper. 
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